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Abstract. Using three-dimensional particle-in-cell simulations, we show that few-MeV 
electrons can be produced by tightly focusing few-cycle radially-polarized laser pulses in a 
low-density atomic gas. In particular, it is observed that for the few-TW laser power needed 
to reach relativistic electron energies, longitudinal attosecond microbunching occurs naturally, 
resulting in femtosecond structures with high-contrast attosecond density modulations. The 
three-dimensional particle-in-cell simulations show that in the relativistic regime the leading 
pulse of these attosecond substructures survives to propagation over extended distances, 
suggesting that it could be delivered to a distant target, with the help of a properly designed 
transport beamline. 
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1. Introduction 

Optical electron acceleration is an important 
application of high-power lasers, in particular for 
the development of new electron and radiation 
sources for probing ultrafast structural dynamics 
in matter [1,2]. So far, research has mainly 
focused on the lower and upper limits of the 
energy spectrum: subrelativistic energies for 
ultrafast electron diffraction (UED) [3-9] and 
ultrarelativistic energies for the development 
of compact radiation sources [10-13]. In 
comparison, little attention has been paid to 
optical electron acceleration in the few-MeV 
range. 

MeV electrons are of interest for relativistic 
time-resolved electron diffraction [14, 15], as 
well as for controlled injection in staged laser 
Wakefield accelerators [16, 17] and dielectric 
acceleration structures [7, 18, 19] that in turn can 
drive table-top synchrotron radiation sources [10- 
13]. In particular for diffraction with ultrashort 
electron pulses, the suppression of space-charge 
effects at relativistic energies allows packing 
more electrons per bunch while maintaining 
ultrashort duration at a distant target. This 
enables single-shot studies of ultrafast structural 
processes that are out of reach of conventional 
(non-relativistic) UED [20]. 

There actually exist ultrafast sources of MeV 
electrons. For example, radio frequency (rf) 
photoguns can produce few-MeV pulses with 
100-fs duration [2, 15,21] and there are proposals 
to obtain even shorter durations (~10 fs) with 
an rf compressor [22]. However, using a laser- 
based acceleration mechanism has the following 
potential advantages: (i) the short wavelength 
of the accelerating field may lead to electron 
bunches with initial durations of the order of 10 fs 
or less; (ii) the intrinsic synchronization between 
the laser and electron pulses is essentially jitter 
free; (iii) compared with nanoclusters and thin 


film targets, optimization using a gas medium 
does not require sub-wavelength alignment of 
the target position; (iv) using a gas medium, 
the target is self-regenerating and can thus be 
used for experiments at high repetition rates. 
Breakthrough laser-driven plasma acceleration 
experiments demonstrated the production of fs 
electron pulses in the sub-lOO-MeV range [23,24] 
and a method to obtain 15-fs pulses in the 5-10 
MeV range was proposed [25]. 

In this article, we show how ultrashort MeV 
electron pulses can be produced by tightly focus¬ 
ing few-TW radially-polarized femtosecond laser 
pulses (RPFLPs) in a low-density gas. We in¬ 
vestigate the laser-gas interaction and electron 
acceleration in a configuration similar to the 
experiment reported by Payeur et al. [5] with 
three-dimensional particle-in-cell (3DPIC) simu¬ 
lations. In particular, we observed attosecond mi- 
crobunching within a longer femtosecond struc¬ 
ture for a laser power in the few-TW range. This 
is about two orders of magnitude lower than pre¬ 
viously demonstrated for direct-field electron ac¬ 
celeration [26,27]. We show that some of these 
microbunches survive to propagation in space, 
suggesting that the proposed method could enable 
the delivery of sub-femtosecond electron pulses 
to a distant target with the current laser and beam¬ 
line technology. 

This paper is divided as follows. First in 
Sec. 2, we present our 3DPIC modelling ap¬ 
proach to relativistic RPFLP-gas interaction. The 
analysis includes an exact analytical solution to 
Maxwell’s equations for the spatiotemporal dis¬ 
tribution of tightly-focused RPFLPs and the com¬ 
plete ionization and plasma dynamics starting 
from a neutral atomic gas. In Sec. 3, we ex¬ 
plain why MeV electron energies and attosecond 
bunching are expected in the few-TW regime un¬ 
der tight focusing conditions. Analytical predic¬ 
tions are confirmed with 3DPIC simulations. In 
Sec. 4, we analyze the propagation of the elec- 
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tron pulses and explain how the electron distribu¬ 
tions could be filtered to isolate the leading sub¬ 
structure. In Sec. 5, we discuss the quality of that 
substructure and its distribution in phase-space to 
show that it is appropriate for being temporally 
compressed and transversely focused. Finally, we 
draw conclusions in Sec. 6. 


2. Relativistic laser-gas interaction model 


To simulate the complete ultrafast and ultra¬ 
intense RPFLP-gas interaction, we performed 
3DPIC simulations with full atomic ionization 
and charge dynamics. In particular, we used a 
customized version of the EPOCH PIC code [28], 
where the electromagnetic field of the RPFLP 
is introduced using the FDTD scattered-field 
formulation [29]. This approach allows us to 
provide the most accurate description of tightly- 
focused RPFLPs via the analytical vacuum-field 
solution. 

For the vacuum field representing the inci¬ 
dent RPFLP, we use an exact closed-form spa- 
tiotemporal formulation given by the following 
expressions in cylindrical coordinates, with r = 
(r,d,z) [30,31]: 
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The physical fields E and B are obtained by taking 
the real part of the corresponding expressions. Eq 
is the amplitude of the transverse electric field 
component and c is the speed of light in vacuo. 
The complex variables are defined as follows: 
R = [r^ (z + cos© = {z + jd)IR, and 

Gf = d’lUit-) ± /(G)] with h = t± RIc + 
jalc. The confocal parameter a is associated with 
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Figure 1. Schematic representation of the laser-gas 
interaction volume analyzed with 3D particle-in-cell 
simulations (3DPIC). 



the Rayleigh range zr by IcqZr = + (koaf - 

1 [32], where ko = InjXQ, do being the central 
wavelength of the laser. 

A Poisson-like pulse spectrum [33] was 
assumed, with the corresponding temporal form 
f{t) = e~-’‘^'> (1 - jojQt! where /o is the 
constant pulse phase and oq = ck^ the central 
angular frequency of the pulse spectrum. The 
parameter 5 is a positive constant related to 
the pulse duration. In the limit where 5' » 
1, the usual carrier-envelope temporal function 
f{f) ^ IT^ ^j(ojot-<j)o) js recovered, with T ^ 

V^/coq. From an experimental point of view, 
the fields given by Eqs. (l)-(3) can be produced 
by focusing a collimated radially-polarized input 
beam with a parabolic mirror of large numerical 
aperture [5,34]. 

Typical numerical simulations are done with 
a cartesian mesh of size 600 x 600 x 2100 
with Na = 30 grid points per wavelength, 

spanning the region x,y e [-lOdo. lOdo] and z € 
[-10/lo5 60do]- We have observed that increasing 
the grid resolution to Na = 40 did not change 
the results significantly. At t = 0, the laser 
pulse is outside the simulation domain and the 
gas is represented by neutral atoms randomly 
distributed in the interaction zone and modeled 
by some 10^ particle markers. A schematic 
illustration of the simulation setup is given in 
Fig. 1. 

Throughout the simulations, multiphoton, 
tunnel, and barrier-suppression ionization were 
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taken into account, as well as the motion of the 
ions and its elfect on the electron trajectories 
(implementation details are found in [35]). Due to 
the low density of the ionized medium, collisions 
and electron impact ionization were neglected. 

We recall that in analogy to the standard 
normalized transverse vector potential parame¬ 
ter ao = ejEj/meCtoo [36], when considering an 
RPFLP it is useful to introduce a normalized lon¬ 
gitudinal parameter a^ = e\E^\peak/fneC(jOo [37], 
where \E^\peak is the amplitude of the longitudinal 
electric field component given at Eq. (2). The rel¬ 
ativistic regime associated with » 1 is charac¬ 
terized by sub-cycle electron microbunching, ex¬ 
plained below (see also [26,37]). 

In direct-field acceleration, the direction of 
the initial longitudinal push to the electrons 
depends on the sign of the local electric field. 
Assuming an infinite homogeneous electron gas, 
half of the free electrons will start moving along 
the laser pulse propagation vector, whereas the 
other half will move in the opposite direction. 
Effectively, the electrons in each field cycle 
will converge locally to the optimal acceleration 
phase position. In the weak field limit {a^ <sc 
1), the expected periodic modulation of the 
longitudinal electron density is imperceptible. 
Nevertheless, in the ultra-relativistic regime 
(aj. » 1), the modulation is strong and 

seen as a series of microbunches, separated 
by Aq. When the microbunches are locked to 
the phase of the laser pulse and accelerated to 
relativistic energies, they are further compressed 
by the local longitudinal electric field, leading to 
microbunches of attosecond duration. 

3. Direct-field electron acceleration and 
attosecond bunching with RPFLPs tightly 
focused in a low-density gas 

In previous publications [8,9], it was shown that 
femtosecond few-hundred-keV electron pulses 


with percent-level energy spread can be produced 
by tightly focusing few-mJ RPFLPs in an hydro¬ 
gen gas whose density is in the 10^^ m“^ range. In 
particular, it was demonstrated that the laser pulse 
parameters and gas density can be optimized to 
cover the 100 - 300 keV energy window that char¬ 
acterizes ultrafast electron diffraction imaging ex¬ 
periments [20]. 

At a dominant wavelength of do = 800 nm, 
typical pulse parameters were koa = 20, s = 70, 
and P = 300 GW, corresponding to a few-mJ 
8-fs full width at half maximum (FWHM) pulse 
focused to a 2Ao spot diameter. These parameters, 
defined by an extensive set of optimization 
calculations [9], were found to offer a good 
compromise between transverse and longitudinal 
effects, with a positive impact on electron energy, 
bunch charge, divergence, and energy spread. 

In the regime explored in [8, 9], the 
normalized longitudinal parameter a^ is slightly 
larger than unity. In such conditions, the 
longitudinal electric field is too weak to induce 
a strong longitudinal compression and attosecond 
bunching. Instead, the ultrashort electron pulse 
duration is due to the acceleration of a thin disk 
of electrons located in a very restricted region 
of the infinite gas target [8]. Electrons outside 
this thin disk region are either deflected away 
from the optical axis or gain little energy from 
the laser field. This effectively results in electron 
pulses with initial sub-femtosecond duration. 
However, the duration of these sub-relativistic 
pulses increases rapidly as they propagate due to 
electrostatic repulsion and the finite initial energy 
spread. 

To predict the conditions to reach the rel¬ 
ativistic regime—where longitudinal compres¬ 
sion would become important—we fitted results 
from [9] for the electron energy as a function of 
the power P of the incident laser pulse. We used a 
^^P scaling law that characterizes direct-field ac¬ 
celeration (see, e.g., [38]). Results are shown in 
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Figure 3. 3DPIC analysis of the transition to the relativistic regime of direct-field electron acceleration by 8-fs (FWHM) 
RPFLPs tightly focused in a low density helium gas (koa = 20, s = 70, 0o = and uhq = 10^^ m“^). (a)-(e) are 2D cuts 
of the 3D electron density structures away from the RPFLP focus along the z coordinate for different values of the laser 
power. (f)-(j) are ID cuts at v = 0 of the corresponding top graphs. Looking at (a) to (c) and (f) to (h), one sees clearly 
the transition from the subrelativistic slicing regime to the relativistic attosecond bunching regime, where the density 
distribution exhibits a sharp modulation at the leading edge due to sub-cycle compression of the accelerated electrons. 
Below the relativistic compression threshold in (a)-(b), the electron density has a peak value less than the initial bound- 
electron density (rie = 2^He)- However in (c)-(e) above threshold, the compressed electron densities have peak values 
greater than 2^He- Still in (c)-(e), it is observed that the peak density increases with the laser power, without important 
changes in the macrobunches structure. Ultimately in (e), the peak charge density exceeds that of (a) by an order of 
magnitude. Total charge in the 3D macrobunches was found to be approximately (a) 1 fC, (b) 2 fC, (c) 3 fC, (d) 4 fC, and 
(e) 5 fC. Charge in the leading bunch is about (b) 0.35 fC, (c) 0.67 fC, (d) 0.85 fC, and (e) 1.05 fC. 



Laser power P (TW) 


Figure 2. Analysis of electron acceleration in tightly- 
focused RPFLP as a function of the laser power. Red dots 
(•) are 3DPIC data points from Fig. 5(a) of [9], while the 
solid black line (-) is a fit of the form a -\- b ^^P (see also 
[38]). We identified in gray two different ultrashort electron 
pulse regimes: the subrelativistic slicing regime, where 
only a thin disk of electron is accelerated (see, e.g., [8]), and 
the relativistic bunching regime, where electrons are axially 
compressed into attosecond pulses (see, e.g., [37]). 


Fig. 2. Conclusions are: (1) average MeV en¬ 
ergies should be reached when P > 2TW and 
(2) for strong focusing conditions {zr ~ Aq), in 
the few-TW power range, the longitudinal elec¬ 
tric field has a super-relativistic strength 
1). According to the ID longitudinal direct-field 
acceleration model [37], longitudinal attosecond 
bunching should then occur naturally. 

To validate our predictions, we simulated 
electron acceleration in helium with the 3DPIC 
method discussed in Sec. 2. For direct compari¬ 
son with the results obtained in the subrelativistic 
electron regime [8,9], we used the following pulse 
parameters: koa = 20, s = 70, and with a 

laser power in the 1-5 TW range. 

Results from the 3DPIC analysis are shown 
in Fig. 3. They reveal that a high-contrast 
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Figure 4. The broad energy spectra associated with the 
femtosecond electron structures shown in Fig. 3. 


attosecond modulation of a longer femtosecond 
electron pulse structure is effectively produced 
when the laser power is in the few TW range. This 
corresponds to sub-lOO-mJ laser pulse energy. 
We recall that previous demonstrations assumed 
Joules (few-cycle pulses, 100 TW) or more to 
produce similar attosecond bunching [26, 27]. 
In these studies, nanoclusters were considered 
to restrict the initial volume occupied by the 
free electrons. Here, the attosecond electron 
microbunches are produced in an infinite gas 
target. This configuration is much simpler, as it 
does not require particular target alignment. 

From the spectra given in Fig. 4, it is 
evident that the femtosecond electron structures 
shown in Fig. 3 are composed of electrons with 
a broad distribution of energy. They are thus 
unlikely to remain together over an extended 
period of time. However, the spectral analysis of 
particular substructures (not shown) suggests that 
these complex femtosecond pulses are effectively 
composed of individual subpulses. We analyze 
this aspect in the following section by looking at 


the evolution of the macrobunches over time, e.g., 
during propagation to a distant target. 

4. Analysis of the propagation to a distant 
target 

Electrostatic repulsion and electron momentum 
dispersion cause the spatio-temporal broadening 
of an electron pulse in time (see, e.g., [39]). 
For ultrashort pulses, the effect is predominantly 
longitudinal, causing the increase of the duration 
as the bunch travels in space [40]. In case of 
a broad spectrum, like that shown in Fig. 4, an 
electron pulse might literally break up before it 
can be delivered to a distant target, e.g., where it 
would interact with a sample to produce a useful 
diffraction pattern. 

We analyzed the propagation of the electron 
pulse structures shown in Fig. 3 with 3DPIC by 
using the moving window technique, were the 
numerical mesh is translated at the speed of light 
to follow the electron pulses as they move. The 
maximum physical time and distance that can be 
simulated essentially depend upon the available 
computational resources. With the parameters 
given in Sec. 2, but where x,y e [-S/to^fido] 
to reduce the computational load, simulating a 
700-fs propagation takes about 40 hours on a 
supercomputer using some 128 processes. It 
corresponds to a distance of 0.2 mm away from 
the RPFLP focus, comparable to the longitudinal 
distances reached in similar PIC studies [25]. 

The results given in Fig. 5 show that at 
the relativistic threshold (2 TW), the initial 
density modulations are lost during propagation. 
However, above relativistic threshold (> 2 TW) 
the leading substructure survives. For the 
electrons accelerated by the 5-TW RPFLP, it 
even maintained a sub-femtosecond duration. 
In comparison, a subrelativistic bunch with a 
comparable initial duration and charge would 
have a duration in the tens-of-fs range (see. 
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Figure 5. 3DPIC analysis of the evolution of the relativistic electron macrobunches produced by a few-TW 8-fs (FWHM) 
RPFLP tightly focused in a low density gas (koa = 20, s = 10, and n = 10^^ m“^). Color plots show 2D cuts of 

the 3D electron density distributions ite away from the optical focus, (a)-(d) zo = 45 jum and (e)-(h) zo = 0.2 mm. (i)-(l) are 
ID cuts at V = 0 of the same distributions: 45 jum as a solid black line and 0.2 mm as a gray area. All distributions were 
normalized for direct comparison. For P = 2TW [(a), (e), and (i)], the structure spreads out and the initial longitudinal 
modulation disappears. For P = 3TW [(b), (f), and (j)], the leading attosecond pulse is still apparent after propagation, 
but with a much broader base. For P = 5 TW [(c), (g), and (k)], the trailing structure of slow electrons spreads out, but 
broadening of the leading pulse base is limited. For P = 5 TW, comparison between simulations with helium [(c), (g), and 
(k)] and argon [(d), (h), and (1)] shows that the formation of the leading bunch does not depend on a specific gas element. 
The initial bunch charge in (d) was about 37 fC, 7 times higher than that of the bunch shown in Fig. 3(e). 


e.g., [9]). 

The macrobunches produced in the longitu¬ 
dinal compression regime share these common 
features: a fast and sharp leading structure fol¬ 
lowed by a slower trailing complex. We have 
observed that by applying a high-pass energy fil¬ 
ter, the leading pulse structure is isolated. Con¬ 
ceptually, such an energy filter is implemented 
by deflecting the electrons with a magnetic field 
and putting an obstacle in the path of the elec¬ 
trons with energies that need to be removed. This 
approach was demonstrated in the femtosecond 
regime by filtering the electron energy directly 
into a magnetic compression device [41]. In prin¬ 
ciple, the same technique can be used for attosec¬ 
ond electron pulses (see, e.g., [42]), although the 
efficiency of the approach in the multi-electron 
regime still needs to be demonstrated. 


To simulate the impact of an energy filter, we 
applied a numerical energy cutoff to the bunch 
structure shown in Fig. 5(h), by rejecting all the 
electrons with an energy less than 5 MeV. The 
resulting bunch, shown in Fig. 6, has a 480-as 
FWHM duration, 0.2 mm away from the optical 
focus. This unexpectedly short duration cannot 
be explained by relativistic effects alone. 

Using the average energy and energy spread 
extracted from the 3D electron bunch of Fig. 6, 
respectively E = 10.2 MeV and dE = 

4.3 MeV, we calculated the pulse duration 
predicted by the relativistic single-electron model 
presented in appendix. With L = 155 pm 
[the distance between the two frames shown in 
Figs. 5(d) and 5(h)], it gives a pulse duration 
of approximately 1 fs, which overestimates the 
3DPIC prediction significantly. We conclude that 
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Figure 6. Spatio-temporal distribution of the macrobunch 
of Fig. 5(h) after electrons with energy less than 5 MeV 
where filtered out. A 2D cut of the 3D electron density 
structure as well as ID longitudinal and transverse cuts 
are shown in (a), (b), and (c), respectively. In (b) it is 
shown that even if the bunch is effectively 0.2 mm away 
from the optical focus, the FWHM duration is still in the 
attosecond range. Calculated bunch charge is 1.5 fC and 
average energy is 10 MeV. 



z!^ 



xIAq 


Figure 7. Normalized (a) longitudinal and (b) transverse 
phase space distribution associated with the filtered electron 
pulse shown in Fig. 6. Both the longitudinal and transverse 
distributions exhibit a nearly linear chirp that is appropriate 
for pulse compression and focusing. 


the space-charge fields from the highly charged 
trailing complex, which are neglected in the 
single-electron model, play a beneficial role in the 
current context by constraining the expansion of 
the leading attosecond substructure at the trailing 
edge. 

Simulating macrobunch propagation with 
3DPIC up to 10 cm, a typical distance in 
UED experiments [20], is impractical (~100 days 
with 1000 computer processes). However, the 
relativistic single-electron model for the bunch in 
Fig. 6 gives an upper limit estimation of 350 fs. 
In comparison, a few-hundreds-keV pulse with 
a percent-level energy spread would be one to 
two orders of magnitude longer (see Eq. (A.2) 
and, e.g.. Fig. 13 in [9]). Still, the phase- 
space analysis we present next suggests that the 
energy distribution in the bunch is appropriate 
for temporal recompression, with the potential of 
sub-fs duration at realistic target distances. 


5. Analysis of the attosecond pulse quality 

We have seen in the previous section that the lead¬ 
ing attosecond pulse shown in Fig. 6 has a large 
relative energy spread {dEjE ^ 0.42). How¬ 
ever, a comparison between Eqs. (A.l) and (A.2) 
indicates that the temporal spreading associated 
with this relativistic pulse is comparable to that 
of a 100-keV pulse with dKjK « 10“^, which 
is very difficult to achieve {K stands for the ki¬ 
netic energy). Close inspection of the longitu¬ 
dinal phase-space distribution given in Fig. 7(a) 
moreover shows that the energy is distributed al¬ 
most linearly in the bunch, which is ideal for re¬ 
compression. The distribution actually suggests 
that the duration of the leading structure after re¬ 
compression could even be shorter than 480 as. 

The thermal transverse emittance associ¬ 
ated with the transverse phase-space distribu¬ 
tion shown in Fig. 7(b) is = 1.5 x 

10“^ mm mrad, similar to known values for ra- 
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dio frequency photoinjectors (typically 6n,x < 5 x 
10“^ mm mrad) [43]. However, due to the large 
angular spread imprinted on the electron distribu¬ 
tion by the curved phase front of the laser pulse, 
the transverse coherence length is small Lj = 
Tij sjipl) ~ 10“^ nm). For a demonstration UED 
experiment, the electron pulse would then need 
to be focused to have a coherence length that is 
sulRciently long to form useful images; typically 
Lt > 1 nm is needed [20]. We emphasize that 
focusing is necessary to keep a maximum of the 
charges in the bunch. But this operation is likely 
to increase the bunch duration, as the net focusing 
elfect of the current magnetic lenses used in UED 
beamlines depends on the electron energy. With 
focusing in mind, it is probably safer to say that 
the final electron pulse duration would be in the 
few-fs range. 

6. Conclusion 

Using three-dimensional particle-in-cell simula¬ 
tions, we have shown that relativistic electron 
pulses can be produced by tightly focusing few- 
cycle radially-polarized laser pulses in a low den¬ 
sity atomic gas. In particular, we observed that 
for a laser power in the few-TW range, longitu¬ 
dinal attosecond microbunching occurs naturally, 
resulting in femtosecond structures with high- 
contrast attosecond density modulations. 

It was demonstrated that in certain condi¬ 
tions the leading attosecond substructure survives 
to propagation over millimeter-scale distances. 
Combined with a high-pass energy filter, focus¬ 
ing element, and dispersion compensator, the pro¬ 
posed method could allow relativistic electron 
diffraction with a temporal resolution in the few- 
fs range, potentially sub-fs. 

Optimization of the laser pulse and gas pa¬ 
rameters is currently being studied to improve the 
final bunch quality and total charge. Preferably, 
future work should include realistic modeling of 


the beamline to include the effects associated with 
the energy filter and focusing element. We em¬ 
phasize that increasing the bunch charge is cur¬ 
rently the most critical issue toward real applica¬ 
tion of the method. Stability of the bunch proper¬ 
ties with respect to shot-to-shot fluctuations of the 
pulse energy and carrier-envelope phase should 
also be addressed. 

Acknowledgments 

This research was supported by the Natural 
Sciences and Engineering Research Council of 
Canada (NSERC). Computations were done on 
the supercomputer Briaree from Universite de 
Montreal, managed by Calcul Quebec and Com¬ 
pute Canada. The operation of this supercom¬ 
puter is funded by the Canada Foundation for 
Innovation (CFI), NanoQuebec, RMGA and the 
Fonds de recherche du Quebec - Nature et tech¬ 
nologies (FRQ-NT). The authors gratefully ac¬ 
knowledge the EPOCH development team. The 
EPOCH code used in this research was devel¬ 
oped under UK Engineering and Physics Sci¬ 
ences Research Council grants EP/G054940/1, 
EP/G055165/1 andEP/G056803/l. 

The authors thank the reviewers for their 
insightful comments and suggestions as well as 
Pran 9 ois Legare, Jean-Claude Kieffer, and Martin 
Centurion for stimulating discussions. 

Appendix 

An electron with velocity v at z = 0 will arrive 
at a remote position z = L at time t = L/v. 
Using the usual relativistic formulas v = (fipjE 
and E = ^c^p^ + m^(fi, this time can be written 
as t = (L/c) x{E/ - m^c^). Then, the arrival 
time delay associated with a collection of non¬ 
interacting electrons whose energies are within 
a certain energy range dE around E is, in the 
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relativistic limit (E » mcfi: 
L 


dt = 


2 4 

me 


dE^^^dE. 
c E^ 


(A.1) 


c [ (E^ - _ 

An estimation of the pulse duration at the target 
due to the finite energy spread is then At + dt, 
where At is the initial pulse duration. It is 
straightforward to show that the non-relativistic 
form of Eq. (A.l) {K <sc mc^) is: 


dt= 

where K stands for the kinetic energy. 
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